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A goal of these studies was development of the measurement methods and algorithms 
necessary to detect wake vortex hazards in real time from either an aircraft or ground-based 
hyperspectral Fourier Transform Spectrometer (FTS). This paper provides an update on 
research to model FTS detection of wake vortices. The Terminal Area Simulation System 
(TASS) was used to generate wake vortex fields of 3-D winds, temperature, and absolute 
humidity. These fields were input to the Line by Line Radiative Transfer Model (LBLRTM), 
a hyperspectral radiance model in the infrared, employed for the FTS numerical modeling. 
An initial set of cases has been analyzed to identify a wake vortex IR signature and signature 
sensitivities to various state variables. Results from the numerical modeling case studies 
will be presented. Preliminary results indicated that an imaging IR instrument sensitive 
to six narrow bands within the 670 to 3150 cm -1 spectral region would be sufficient for 
wake vortex detection. Noise floor estimates for a recommended instrument are a current 
research topic. 


I. Introduction 

T He overall objective of this work was to investigate the feasibility of detecting several types of in-flight 
aviation hazards with a Forward Looking Interferometer (FLI): clear air turbulence (CAT), wake vortices, 
dry wind shear, icing, hazardous runway conditions. 1 The FLI is a 320 x 256 pixel imaging FTS that operates 
in the IR band (650 cm~ l to 3300 cm~ l wavenumber). 2 The FLI concept combines high spatial, spectral, 
and radiometric resolution and leveraged remote sensing techniques developed for satellites. 3 The FLI is an 
instrument concept that, in the long-term, would obtain in-flight measurements that would alert flight crews 
to potential weather hazards during all phases of flight. 4 The research described here was funded by the 
NASA Aviation Safety program, Atmospheric Environmental Safety Technologies project. 

A goal of this study was the development of the measurement methods and algorithms necessary to detect 
wake vortex hazards in real time from either an aircraft or a ground mounted FLI instrument. The Terminal 
Area Simulation System (TASS) was used to generate wake vortex fields of 3-D winds, temperature, and 
absolute humidity. 5 These fields were input to the Line by Line Radiative Transfer Model (LBLRTM), a 
hyperspectral radiance model in the infrared, employed for the FTS numerical modeling. 6 

Extensive numerical modeling was conducted using a combination of the two models, TASS and LBLRTM. 
An initial set of cases was performed to identify wake vortex IR signature sensitivities to various state 
variables. An algorithm was developed to analyze the three-dimensional data sets that were generated by 
the combined models. This algorithm selects the FTS channel with the maximum signal, and enables data 
reduction to two-dimensional imagery. Results from the algorithm are used to determine a wake vortex 
thermal signature and to simulate FTS imagery. 
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II. TASS Model Parameters 


For this research, the TASS data sets consisted of three-dimensional grids extracted every 10 seconds 
after vortex insertion and continuing for a total duration of 60 seconds. Table 1 lists the relevant parameters 
used by the TASS model. 


Table 1. TASS model initialization parameters listed 
below were used to generate Boeing B74T wake vor- 
tices at 10 second intervals after initialization. 


Parameter Boeing B747 

Grid (X by Y by Z) 800 m by 600 m by 600 m 

Initial vortex separation 50 m 


TASS model circulation was determined by the 
type of aircraft selected and assumed an aircraft with 
a clean configuration in level flight. Additional pa- 
rameters are listed in Table 1. In the table, the mod- 
eled grid had dimensions A, or along track; T, or 
cross-track; and Z, or vertical. This spatial nomen- 
clature was used throughout this research. 


Initial vortex height 500 m 

Initial vortex core radius 4.5 m 


III. IR FTS Numerical Model 


For this modeling effort, the LBLRTM was used 
to compute radiances within a specified atmosphere 
of various gases over a specified IR band. 7,8 The 
LBLRTM code required input parameters describing 
atmospheric gas temperatures, pressures, and con- 
centrations. 9 An optical ray, originating at the top 
of the atmosphere extended through the TASS mod- 
elled wake vortex to the IR FTS instrument mounted 
on a following-aircraft. Each FTS pixel was simulated 
as the collected IR energy along a single optical ray. 
Each point along the ray was defined by Z- axis height, zenith angle to local vertical, and temperature and 
moisture values of interpolated points embedded in the U. S. Standard Atmosphere. 10 

A list of LBLRTM parameters and ranges of values is given in Table 2. For these simulations, the 
unapodized spectral resolution was configured as 0.5 cm -1 . The output data was apodized using the Beer 
function 11 to yield output with various spectral resolutions. For the results presented, the wavenumber 
spectrum ranged from 650 cm -1 to 3300 cm -1 at 1 and 16 cm -1 intervals. Simulations at lower spectral 
resolutions were performed but are not presented. 


Initial vortex circulation 

565 m 2 /s 

Initial vortex descent velocity 

1.8 m/s 

Initial temperature (at 0 m) 

288.15 °K 

Initial pressure (at 0 m) 

1013.25 hPa 

Temperature lapse rate 

7.4 °K/Km 

Brunt-Vaisala frequency 

0.0015 

Eddy dissipation rate 

4x 10 -5 m 2 /: 

Water vapor mixing ratio, Q v o 

6 g/kg 


Table 2. All simulations were conducted using the 
LBLRTM simulation parameters listed here. Typi- 
cally, only one of these parameters was varied at a 
time. 


Parameter 


Range or Value 


The altitude of the aircraft and wake vortex was a 
parameter with values ranging from 10 Km (a typical 
cruise altitude), 3 km, a typical holding pattern alti- 
tude, approach or 1 km, to landing or 0 km. Wake 
vortices may be persistent at any altitude and can 
be encountered by other aircraft regardless of flight 
phase (takeoff, cruise, hold, approach). 

The LBLRTM model was configured to use the 
following atmospheric gases: H 2 O, CO 2 , CO, NO 2 , 
NO, CH 4 , and SO 2 . The wake vortex was mod- 
eled as a region along a ray of elevated temperature 
and diminished absolute humidity. This model as- 
sumed that a generated vortex pair captured cooler, 
dryer air from a higher altitude, and descended into 
a warmer, moist er ambient environment. As a worst 
case scenario, the leading and following aircraft were collinear, such that the vortex pair posed the greatest 
hazard. In this case, the FTS field of view has rays aligned with the vortex axis. The vortex was 500 m 
from the following aircraft and had a depth along the ray of 300 m. The elevation angle corresponded to 
FTS instrument pointing angles distinct from aircraft angle of attack. 


Spectral Resolution 
Range of IR Band 
Aircraft Altitudes 
Elevation Angle 
Distance to Vortices 
Wake Vortex Depth 


16 cm -1 (apodized) 

650 cm -1 to 3330 cm -1 
0 km, 1 km, 3 km, 6 km and 10 km 
0°, 5°, 10°, and 15° 

500 m 
300 m 
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IV. Maximum Channel Standard Deviation Algorithm 

The following algorithm was devised to reduce both simulated and measured FTS data, though only the 
simulated data is analyzed in this research. The algorithm performs a statistical comparison of all pixels for 
a particular wavenumber channel, in order to select that channel with the most signal. 

Within each FTS data cube Tb{u , z,v), a single channel contains an array of brightness temperature data 
in the two spatial dimensions, y and z, and spectral dimension, v. Let L and M be the number of pixels in 
the horizontal and vertical dimensions, respectively. Also, let N be the number of channels, or wavenumbers 
in the spectral dimension. The first step in this algorithm eliminates natural radiance variations due to 
elevation angle. Brightness temperature deviations removed the row mean using 

1 M 

T* B (i,j,k) = T B (i,j, k) - (1) 

3 = 1 

for i = [1 . . . L], j = [1 . . . M], and k = [1 ... TV]. The next step was to compute the spatial mean from the 
brightness temperature deviations. The spatial mean was defined as 

L M 

^) = CTEE^ (i -i^) (2) 

i= 1 3 = 1 

The spatial standard deviation for each channel was computed as 


°Yz{k) = 



L M 


i=l j = 1 


(3) 


On a per channel basis, the spatial standard deviation was used to define a thermal signature for a wake 
vortex. Results from simulations that varied altitude are shown in the next section. 

The next step was to search for the channel, Vmax , with the maximum standard deviation. 


Vrnax — k = max [a Y z{k)\ (4) 

k 

where N was the number of channels and k was wavenumber. By this process, the spatial brightness 
temperature at v m ax was selected for each pixel, forming a Maximum Channel Standard Deviation (MCS) 
image. 


MCS(iJ)=T^iJ,u rnax ) (5) 

Results are shown as images of Tb at Vmax f° r simulated FTS imagery. 

V. Simulation Results 

TASS output for the Boeing B747 wake vortex temperature field is shown in Figure 1 (a) and absolute 
humidity field in Figure 1 (b). These images are from 3-D data in the YZ plane at X=0. The thermal 
variation illustrated in the data was about 1 °K, using the color scale for dark red to yellow. This small 
temperature variation is largely a function of the temperature lapse rate of the atmosphere. 

As described in Section III, TASS data was embedded within the Standard Atmosphere. A set of rays 
through this combination were used to simulate FTS pixels. Each ray corresponded to a single pixel. The 
geometry of this scheme was parameterized and could be varied to study various effects. The primary variable 
of interest for this research was altitude. 

The output LBLRTM radiance estimates exhibited natural variation unrelated to atmospheric effects. 
The Planck function was used to remove spectral variability. A typical brightness temperature spectrum 
for a single ray extending from the top of the atmosphere down to the FTS is shown in Figure 2. For this 
particular plot, the FTS was simulated at an altitude of 1 km and with an elevation angle of 5 degrees. The 
range of temperatures is indicative of the emission spectrum and is too great an amplitude to discern the 
small variations associated with a wake vortex. 
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Figure 1. TASS Boeing BT47 at 60 s after initiation showing wake vortex fields. 
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Figure 2. Sample brightness temperature spectrum from LBLRTM at 1 km altitude and elevation 5 degrees. 


The wake vortex thermal signature is illustrated in the following plots of ayz, shown for FTS simulations 
with 16 cm -1 spectral resolution. Figure 3(a) is at altitude 0 km, Figure 3(b) at 3 km, Figure 3(c) at 6 
km, and Figure 3(d) at 10 km. At this lower spectral resolution, the FTS model has smoothed or filtered 
the higher resolution spectral content. The result is a set of plots that reveal the thermal signature of the 
wake vortex. In these plots, the wake vortex thermal signature can be clearly defined as consisting of six 
distinct bands. More importantly, these results illustrate that a high-resolution FTS is not needed, and that 
an imaging radiometer with at least 16 cm -1 can be used instead. 

As shown in Figure 4, the ayz data can be used to select wavenumber bands that represent the wake 
vortex thermal signature. From this plot, the wake vortex thermal signature can be defined as incorporating 
the following five narrow spectral bands: 670 to 750 cm -1 , 1250 to 1450 cm -1 , 1580 to 1660 cm -1 , 1800 to 
2100 cm -1 , and 2200 to 2350 cm -1 . A sixth band above 3000 cm -1 was excluded due to non-monotonic 
altitude variation. From the plots of Figure 3, the different IR bands composing the wake vortex thermal 
signature are seen to exhibit altitude dependence. Bands 1 and 5 are amplitude invariant with altitude, 
while Bands 2, 3 and 4 merge together above 3 km. Band 6 as evidenced by the appearance of a large spike 
in Band 6 of Figure 3(d). By utilizing these bands, an imaging IR instrument will have higher sensitivity 
to the vortex radiative emission and be less sensitive to altitude dependencies exhibited in the plots. The 
specific bands with optimal sensitivity depend upon altitude and atmospheric moisture concentration. 

Results from the MCS algorithm, as defined in Section IV, are presented. A simulated FTS image plot 
with resolution 16 cm -1 at altitude 0 km is shown in Figure 5(a) and at altitude 3 km in Figure 5(b). In 
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Figure 4. Illustration of wake vortex thermal signature with identification of six IR bands. 


these imagery plots, a single channel represent the entire spectrum at each pixel. This channel is selected 
using the MCS algorithm. Note that the wake vortex pair is evidenced by semicircle of warm, moist air. 
Due to long computation times for the large number of rays, subsampling of the original TASS data results 
in the pixelation seen in these figures. The temperature range for these plots is about 0.025 °K using the 
color scale for dark red to yellow. The variation does not degrade significantly with altitude. 

A second set of FTS images is shown in the plots of Figure 6(a) at altitude 6 km and in Figure 6(b) 
at altitude 10 km. As with the previous plots at the lower altitudes, the wake vortices are evident. How- 
ever, these plots illustrate a loss of signal amplitude with increasing altitude. These plots indicate that a 
thermal signature associated with wake vortices is theoretically capable of being detected by an imaging IR 
instrument. 

VI. Conclusion 

The LBLRTM was employed to model FTS responses to atmospheric variations. Several parameterization 
cases were performed to identify wake vortex IR signature sensitivities to various state variables. These 
variables included temperature, water vapor, atmospheric gas, turbine jet exhaust gas concentration, altitude, 
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Figure 5. FTS simulations at resolution 16 cm 1 for IR Band 650 cm 1 to 3350 cm 1 at elevation angle 5°. 



Figure 6. FTS simulations at resolution 16 cm 1 for IR Band 650 cm 1 to 3350 cm 1 at elevation angle 5°. 


viewing angle or elevation, and distance from FTS model to wake vortex. 

The results of this preliminary study revealed that a high resolution FTS-type instrument may not be 
needed to detect a wake vortex thermal spectral signature. Additional simulations are needed to examine 
instrument response to wake vortex distance and resulting imagery generated from the cryz data. 

On a per channel basis, the spatial standard deviation was used to define a thermal signature for a wake 
vortex. The primary variable of interest for this research was altitude. Results from simulations that varied 
distance are a future research topic. For the simulated wake vortices, the wake vortex thermal signature was 
clearly visible in the resulting imagery. 

In these idealized simulations, no instrument noise was modeled. As such, a required noise floor needs 
to be estimated. This necessary step requires the simulation of FTS or radiometer imagery. In addition, 
the spectral response of the IR sensor was not simulated. Typical IR imaging FTS designs employ either 
HgCdTe or InSb focal plane arrays as the IR sensor, each having an associated responsivity spectral band. 
Unfortunately, this research identified five narrow bands that extend across both of these two different IR 
sensor types. 

Based on this research, a recommended instrument type for wake vortex detection is an imaging IR 
radiometer sensitive in the following five narrow spectral bands: 670 to 750 cm -1 , 1250 to 1450 cm -1 , 1580 
to 1660 cm -1 , 1800 to 2100 cm -1 , and 2200 to 2350 cm -1 . Further research is needed to determine a 
recommended a maximum instrument noise floor after additional simulations are performed. After altitude, 
the next major variable of interest is change in distance from the wake vortex. Yet more analysis is required 
to compare these simulated results with measured data. 

Additional research is needed to analyze the impact of turbine jet exhaust on these results. How will 
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rolled-up exhaust gas change the thermal wake vortex signature for these radiative transfer simulations? Can 
the IR chemical spectroscopic signal of turbine jet exhaust be used as a proxy for wake vortex detection? 
The numerical techniques describe here should enable the research to resolve these questions. 
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